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Abstract. 

Scalar mediated interactions among baryons extend well above the Compton 
wavelength, when they are embedded in a Bose-Einstein condensate composed of the 
mediating particles. Indeed, this non-trivial environment results in an infinite-ranged 
interaction. We show that if the Dark Matter of the Universe is composed of such a 
condensate, the imprints of an interaction between baryonic and Dark Matter could 
be manifest as anomalies in the peak structure of the Cosmic Microwave Background. 
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1. Introduction 

In work carried out over past years, we have studied dispersion forces £Q (i.e. van 
der Waals type forces) exerted among macroscopic bulk matter mediated by light (or 
massless) particle exchange, both in vacuo [21 Ej and in a heat bath Known 
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examples are the Casimir-Polder forces j3] among neutral atoms (where the exchanged 
particles are photons) or the Feinberg-Sucher two-neutrino exchange forces |B] among 
baryonic matter. Because real matter sits in the cosmic microwave background and 
in the relict neutrino background, those forces are affected by the photon thermal 
background and the neutrino background, respectively [U Ej. Now, scalar particles 
are fundamental ingredients of the Standard Model of Particle Physics (SM) and 
completions thereof. Known examples are, Higgs bosons, axions, majorons, scalars 
appearing in supersymmetric extensions of the SM, dilatons, radions, etc. Hence, when 
coupled to ordinary matter, their exchange leads also to dispersion-like forces jHj- 

Scalar particles, such as the axion, have been considered in the past as candidates for 
the dark matter in the Universe. While the favorite candidates for the dark matter in the 
Universe are weakly interacting particles with masses in the GeV range (WIMPS), the 
most popular example of this being the lightest supersymmetric particle (in particular, 
the neutralino) 0, it cannot of course be excluded that a scalar (even one not given in 
the list above) is actually the real ingredient of dark matter [TUj . In this respect, it has 
been recently proposed that the observed 511 keV emission from the Galactic bulge 
could be the product of very light annihilating scalar dark matter particles |12j . 

Therefore, the putative forces caused by the double exchange of such scalars would 
be also affected by the presence of the scalar dark matter. When two scalar particles are 
exchanged in the t-channel, spin-independent dispersion forces are being generated that 
add coherently over unpolarized bulk matter and extend over distances on the order of 
the Compton wavelength of the mediating scalars. If the fundamental coupling of these 
scalars to matter fermions is of the usual Yukawa type, then the long distance behavior 
of the associated potential in vacuum is specifically of the form ~ exp(—2mr)/r 5 ^ 2 |13j . 
When the bodies subject to those forces are embedded in a heat bath made of the same 
scalars, a very dramatic effect takes place provided the scalar particles carry a conserved 
charge and the boson gas reservoir is characterized by a nonzero chemical potential. 
Namely, the finite range potential (now behaving as ~ exp(—2mr)/r 2 ) becomes infinitely 
ranged when the heat reservoir suffers Bose-Einstein condensation. Specifically, an 
~ exp(—2mr)/r 2 potential turns into ~ 1/r. Conversely, when the heat bath makes the 
transition to the uncondensed phase, this infinite range force becomes finite ranged. The 
phenomenon, first described in ref f3j (see also J3]), comes about as a combination of 
kinematics (three-momentum exchange of the matter system with the medium) and the 
collective effect of condensation of charge. In this paper we wish to explore a physical 
realization of this effect. Our purpose is to display a system where the phenomenon 
is of physical import. Since our primary goal is to investigate the implementation and 
physical consequences of Bose-Einstein condensation of a relativistic scalar gas itself, we 
will not focus on any specific particle physics model (see e.g. [T3] for the effects on the 
Standard Model particle dispersion relations) nor claim that the envisaged cosmological 
arena is a realistic one either. We consider it as a convenient playground for testing the 
potentialities of a model characterized by matter in a light scalar heat bath. 

On the other hand, one should be open minded as to which is the actual nature 
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of dark matter and consider other options as long as WIMPS are not experimentally 
detected (see e.g. ^H] for recent accounts of the status of particle dark matter). Following 
the progress of experiments in measuring the cosmological parameters, the so-called A 
CDM scenario has emerged as the Standard cosmological scenario f7j suggesting that 
most of the energy density of the Universe is stored in a dark sector, 2/3 of which, 
dubbed dark energy, has negative pressure and the rest is composed by pressureless 
dark matter (although see [TH] for alternative fits to the experimental data). The 
simplest alternatives are a cosmological constant, A, to explain the recent period of 
accelerated expansion and WIMPS for the dark matter, but other candidates have been 
put forward that usually entail a richer structure in the dark sector. To name just a few, 
dynamically evolving scalar fields are at the basis of quintessence models for the dark 
energy , condensates may play a significant role in many cosmological phenomena [20] 
and provide a unified description of dark energy and dark matter [25, also long range 
forces in the dark matter sector motivated by string theory have been examined in |22j . 

In section [2] we introduce the Bose-Einstein condensate of scalar dark matter 
particles and analytically study their influence on the peak structure of the CMB. A 
numerical example is next discussed in section El before ending with our conclusions in 
section |H 

2. Bose-Einstein condensate of Dark Matter 

2.1. The set-up 

Let us start by displaying the model. We assume the existence of light scalars that 
carry a conserved quantum number Q. These scalars constitute the bulk of dark matter. 
For a relativistic boson gas [221 EH E3 i-e. m/T « 1, condensation occurs below 
T c = (3q/m)^ for a fixed charge density q. Conversely, for given T, when the charge 
density exceeds q c = mT 2 /3, Bose-Einstein condensation follows. Now, in an expanding 
Universe number densities vary with T 3 and hence there should be a period in the 
history of the early Universe where a condensed phase of scalars coexists with a gaseous 
scalar phase: for a sufficiently high temperature, the charge density exceeds q c [23]. The 
scalars, while they are in a condensed phase, are a pressureless fluid and thus constitute 
cold dark matter. In order that structure formation can proceed as usual, we certainly 
want the scalar background to be in the condensed phase at matter-radiation equality 
and beyond, and the eventual transition to an uncondensed phase to take place not 
before the large scale structures we observe today have had enough time to form. We 
shall see below how these conditions can be met. 

We suppose next that our scalars couple to ordinary baryons via the effective 
interaction 



where mj is the nucleon mass, arbitrarily chosen to make the effective coupling g e ff 
dimensionless. It may be viewed as the low energy limit of the interaction shown in 



9eff 




(1) 
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Figure Q involving the fundamental underlying Yukawa interaction, with strength A, of 
the scalars to heavy fermions associated to a high energy scale M. The dimensionless 
coupling g e ff can be understood as g e ff = X 2 (mb/M). Equation (0) leads, via the 
exchange of two scalars depicted in Figure El to the asymptotic potential (i.e. r ^> mT 1 ) 



V 



-9eff T 



647r 2 m^r 2 " ^ ^ 

between two static matter fermions sitting in a scalar heat bath at temperature T > T c , 
where T c is the Bose-Einstein condensation temperature 




1> 




Figure 1. Effective coupling of the scalars <p± to ordinary baryons ip. 




Figure 2. Exchange of two scalars giving rise to the force among baryons. 
Below T c , i.e. in the condensed phase, the resulting potential reads |T3] . 

v ^ _jhf^_ = _j!ifll (3) 

The effective interaction (P) is the simplest form describing the exchange of two 
scalars among baryons. Other forms, like e. g. derivative couplings appropriate for 
Goldstone bosons, would result in a different r behaviour [3] of the potentials (0) and (jSJ), 
but the transition to a long range interaction below T c , with the removal of the 
exponential damping factor, remains unaffected. 

Having mentioned Goldstone bosons, let us, at this point, note that there is 
an interesting parallel between spontaneous symmetry breaking and Bose-Einstein 
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condensation: both types of systems possess bosonic fields which acquire a constant 
classical component. This parallel is not exact, however, since the physics of how these 
constant classical components arise is different (see e.g. ref [23] )• For instance, the 
noninteracting charged scalar field cannot exhibit spontaneous symmetry breaking, but 
it does show the phenomenon of Bose-Einstein condensation and this is indeed the 
process that occurs in our setting. The scalars that are exchanged in Figure El are not 
Goldstone bosons (note e. g. that they have a non-zero mass). Linear couplings are the 
simplest alternative to describe their interaction with the baryons, and thus our choice 
in eq. JT]). 



2.2. Influence on the CMB 

Because of this long range dispersion force between baryons, the Euler equation for 
baryonic matter in the tightly coupled photon-baryon fluid gets an extra body force 
term kS, (in Fourier space; k is the wave-number associated to a sub-horizon sized scale) 
where the potential £ obeys the Poisson equation: 

k 2 ^ = -47TK{a/a ) 2 p b 5 b (4) 

with k = ^f^fr and where a is the scale factor, the subscript refers to "today", and 
Pb and 5b are the baryon energy density and density contrast, respectively. 

Combined with the continuity equation for the photon temperature fluctuation 
= one obtains the oscillator equation responsible for the acoustic peaks in the 
power spectrum of the CMB: 

e + — — e + k 2 c 2 s e = -— * - — — $ - $ (5) 

l + R s 3 l + R v ' 

where \& and $ are the gravitational potentials associated to dark matter and R = |^-. 
The derivatives are with respect to conformal time rj. The key difference with the usual 
acoustic equation [2H1I2ZI is the "speed of sound" . Indeed, instead of the proper velocity 
of sound c 2 = |y^r, we have c 2 = c 2 {l — e s {k)) where, 

127TK 

T 2 

As in the usual case, the equation can be solved very accurately with the help of the 
WKB approximation. The phase of the oscillations is given by ks with 



e s (k) = 3 x IQ^bhY^^ieV/cm 3 ) (6) 



c s dr] (7) 

i.e. the effective sound horizon. 

Modes that verify k n = nn/s rec at recombination are caught at maxima 
(compression) or minima (rarefaction) of their oscillation and render the acoustic peaks 
in the power spectrum. Since multipole moments / oc k approximately, 



5L 5k„ 5s 
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where 5 s = s — s. As a consequence of the infinite-range interaction, the Z-position of 
the peaks will be shifted (except the first one, see below). Treating e s {k) in equation ((7j) 
small perturbation we find: 

5s 



s 



2s(r] r 



e s (k)c s dr] (9) 

Vent 



Note that the above integral starts at r] ent , which is the time at which the k n -mode 
enters the horizon. From this moment onwards, the causal interaction mediated 
by 20-exchange is operative. The integrand depends on time because the speed of 
sound changes with time but also because the parameter e s {k) depends on the critical 
temperature T c which decreases as the Universe expands. The integral is easily done 
analytically but the explicit expression is cumbersome. If we work it out a little bit by 
introducing its numerical value at recombination, we obtain: 



5s 



4.31 x 10 13 #, 



2 
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(k/OMMpc- 
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(0.75 - 7.7AR ent ) -2 



2 

ent 



18.3 In- 



D2 

O.leVV fVthh 2 \ [VL DM h 2 



J V0.02 J \ 0.14 



(10) 



vo.i7 +R en t + o.4iyrr^ cut 

In equation (fTUj) we used the value R rec = 0.46 and for R ent we shall use 

R ent = 3.66 {k/OMMpc-y 2 (^y^) (11) 

which is adequate for scales that enter the horizon after matter-radiation equality. 

Equation (|10|) evaluated at the first peak (i.e. k ~ O.OSMpc" 1 ) vanishes. It is rightly 
so, because the first peak corresponds to a scale that entered the horizon just before 
recombination and therefore the dispersion force has had no time to operate. Higher 
harmonics, that correspond to smaller scales and hence had more time to oscillate, 
give positive shifts (towards larger Fs) because our force effectively reduces the sound 
horizon. Finally and obviously, the higher the multipole, the larger the shift. 



3. A numerical case study 

The height and location of the acoustic peaks depends on the parameters of the working 
cosmological model. For instance, the location of the first peak strongly favors a flat 
Universe with a less than critical baryon density and, combined with the existence of 
the second peak, one has an estimate of the total cold dark matter density [23 I2Z|- 
Modifications of the cosmological model will show up, thus, in the peak structure, 
and constraints on such variations can be obtained by means of numerical likelihood 
analysis. For example, the speed of sound of the dark energy is c s -C 1 in models with 
noncanonical kinetic energy density, such as k essence, and densities of the order of 1% 
can be distinguished from models whose sound speed is near the speed of light |28j . 
Theories of neutrino mass generation can lead to a uniform shift of the peaks to larger 
/ and an enhanced damping due to a delayed matter-radiation equality |2"H] . 
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As we do not want to endeavor a full analysis of the allowed parameter space in our 
model, for the sake of definiteness, we take the uncertainty in the second peak position 
to be ^ = 1%. This figure provides an order of magnitude estimate of present and 
planned experiments' ability to reconstruct the CMB power spectrum |27j . 

Then the third peak moves by a 3%. From equation (|1()J) we infer that these 
variations follow for an arbitrarily chosen scalar mass = O.leV^and for a coupling 
g e ff ~ 2.3 x 10~ 8 , which is a very small number. Let us now check the self-consistency 
of the picture when we adopt this value for g e ff- 

The scalar boson gas is made of particles with +1 Q-charge, which we denote by 
</>+, and anti-particles with —1 Q-charge denoted by </>_. To easy notation, however, 
we shall display the charge sign only when strictly necessary. Let the net total charge 
be positive (this is purely conventional) and therefore there is an excess of particles 0+ 
over antiparticles 0_ in the Universe. At the very early epochs of the Universe, say 
at T > M, when the degrees of freedom associated to the high energy scale M (e.g. 
fermions x with charge Q in Figure are fully relativistic, our scalars are in thermal 
contact with all relativistic species, the density is very high, and the gas of scalars is 
condensed with most of the charge sitting in the condensate. The condensed phase 
consists overwhelmingly of particles 0+ of zero momentum. They are the bulk of dark 
matter. Since charge Q in any co-moving volume is conserved, the charge density q is 
bound to be: 

a, fVL DM h 2 \ f0.1eV\ o , , 

,= 1.5 x 10*(1 + ;)*(^)(— )*n-> (12) 

Although charge in excited states is a smallish fraction 0(m^/T) of the charge in the 
condensed phase, the number density of scalars and anti-scalars occupying excited states 
is comparable to the other relativistic species at those large temperatures. Later on, 
when temperature has dropped below the energy scale M, and consequently Q-fermions 
are extinct, our scalars might still interact with baryons through elastic 0iV-scattering 
(N stands for nucleon) where the 0- interaction rate T ~ rib(<Jv) contains the properly 
averaged cross section x relative velocity: 

M = ^^oV / TT^TTWToM 2 ! 1 + f BE )W5\P - P>) (13) 
2(p2> 2(pg) J (27r) 3 2p u (27r) 3 2p; u 

In this expression \M\ 2 is the spin averaged transition matrix element squared and 
we have included the spontaneous plus stimulated emission factor (1 + /be), where 
/re = [exp (p® — fjLj /T — 1] _1 is the Bose-Einstein distribution function (/x being the 
chemical potential). This factor weighs the population of final particle phase space 
in a boson medium. The calculation of the interaction rate for thermal scalars (i.e. 
those with nonzero momentum) gives a value, at T ~ lGeV, comparable to the Hubble 
expansion rate but drops precipitously beyond iVW annihilation due to the massive loss 
of scatterers. The calculation of the interaction rate for the zero momentum scalars in 
the condensate is not as straightforward. 
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Indeed, the final scalar phase space integral in the previous equation does not 
properly account for the sum over states in the condensed regime. The zero modes are 
not included in the sum. This is not a problem in the non-degenerate case because charge 
is thinly distributed over the states and the ground state gets only an infinitesimal share 
of the total charge. However, in a condensed Bose-Einstein gas a macroscopic fraction of 
the charge sits in the ground state. In our case, the ground state is populated essentially 
by particles <f> + . Actually, the integral (evaluated at fixed chemical potential \i = m^) 
gives then the contribution of the excited states to the phase space. We shall estimate 
the stimulated emission contribution to the interaction rate to be: 

n 

r s ti m ^n b {a*v)— (14) 
n 

where (<r*v) is explicitly given byf: 

where we have omitted the spontaneous emission term which in this case is negligible 
compared to the stimulated emission term, and the fraction n/n* of the total number 
of particles over the particles in excited states can be approximated by: 

2~ 

q 



n n k=0 + n _ ; q k=0 _ ^ 



(16) 



n* n* n* ((3)T 3 /tt 2 

In equation (|15|) we took into account that the rates are to be evaluated in a co-moving 
frame at rest with respect to the Hubble flux. In this frame the scalars in the condensate 
are motionless and the baryons are non-relativistic in the relevant temperature span 
(T ~ lGeV and below). The energy-momentum deltas can be used to crack down the 
phase space integrals and to obtain a much simpler expression that reads: 

= TT^-TT 2 - / - W //n 1 — - P 4 r- (17) 



2m <t> 2m b J 4n exp(p£ - m^)/T - 1 p /o + p io f± + Pb z \ 

with p b — |p 6 |, p' b — \p' b \ and z is the cosine of the "laboratory" scattering angle. The 
kinematics of the scattering process corresponds to an extremely glancing collision as 
required by energy-momentum conservation and the delta function in equation (|17|) 
describes approximately the fact that 1 > z 2 > 1 — (m^/mj) 2 . In equation (|T7|) . p' b and 
p'2 are fixed by energy-momentum conservation: 



and 



m b 



^(1+^4) 



2 m 2 b 

| the star in a* is a reminder that we are summing only over excited final states 
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Collecting the different pieces and introducing them in equation ()17|). we obtain: 

i6^ m ^ UJ 



^) = |4<S 1/2 (20) 



where we used ~ ^fim^T as dictated by the equipartition theorem. 

We get the final form for the interaction rate of scalars in the condensate 
by plugging the equation above into equation ()14|) where the factor n/n* is about 

70 ) i~^f) ir^tm) in this temperature regime (T ~ lGeV and below). Even 
below anti-baryon extinction, when the baryon number density rib is 1-22 x 10 _10 T 3 , the 
resulting rate is many orders of magnitude larger than the Hubble rate: 

I « 4 x io- ( ) 2 ( 3 1/2 (2™>r\ 1/2 3/2 (21) 

H V 2 - 3x10 " 8 / V m * / V m b J V 0.14 J \lGeV J V 7 

Hence, while thermal scalars are already decoupled at temperatures below or about 
lGeV, the zero modes in the scalar condensate are in full thermal contact with the 
cosmic plasma. Though decoupled, thermal scalars keep their equilibrium distributions 
with momenta redshifting as the Universe expands. And hence they keep a common 
temperature with the condensate bosons that are tightly coupled to baryons and 
photons. A common temperature is maintained only until anti-baryon extinction takes 
place. By then photons are reheated with respect to thermal scalars and the condensate, 
still coupled to the baryon-photon fluid, shares its temperature. Actually the zero modes 
in the condensate keep in thermal contact with the baryons well beyond matter-radiation 
equality (see below). A further widening of the temperature gap between thermal scalars 
on the one hand and photons and the condensate on the other hand occurs at two more 
stages: after and later after e + e~ annihilation. Hence at nucleosynthesis the ratio 

of temperatures is: 

T*m.= ( rel.d.o.f.T<m e \ 1/3 

T \rel.d.o.f. T > lGeV J V ; 

Since the effective number of degrees of freedom associated to the thermal scalars at 
nucleosynthesis is proportional to this fraction to the fourth power, there is clearly no 
conflict with primordial helium yields [SD]- We turn next to the question as to when 
does the condensate loose thermal contact with the baryon-photon plasma. 
In a matter dominated Universe the expansion rate is given by: 



^Vir"" (23) 

Pm being the non-relativistic matter density. 

Decoupling of the particles in the boson condensate from the baryon-photon fluid 
happens when their interaction rate roughly equals the cosmic expansion rate. Actual 
numerical comparison of both rates renders: 

J_ ^ / g eff \ 2 (O.leVV (IGeVV 12 ( Sl DM h* \ ( 0.15 V' 2 f T V 

H V 2 -3xlO"V V m 4> ) V m b J V 0- 14 J \tt M h 2 J V0-016eUy 1 } 
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This is a rough estimate that determines the moment of condensate decoupling to 
happen about a factor 10 — 20 below the recombination temperature. Here we took 
into account that scalars in the condensate and thermal scalars do not share the same 
temperature at these later stages. Indeed, in this temperature regime we replaced T^^o) 

for T in equation (|HJ), we used n/n* ~ 70 (^f) (^f) (r^^f ' and we related both 
T and T^^q) via equation (J2*2*j) . 

The charge density at decoupling is still about a factor of five above critical. So 
the transition to an uncondensed phase is yet to happen. However, since scalars are 
non-interacting and thermally isolated from ordinary matter ever since decoupling took 
place, scalars will maintain the equilibrium distributions they had at decoupling. The 
zero modes, i.e. the condensate, will therefore keep in the lowest energy state beyond 
the critical point because no interaction can pump them out of the ground state. As a 
consequence, baryons in the cosmic scalar background today, i.e. at zero redshift, should 
feel the long-range two body exchange forces under scrutiny. They are extremely feeble. 
In fact, being proportional to the present cosmic Q-charge density and to ~ 5 x 10~ 16 
(see Equation (j3J)), they turn out to be almost a million times weaker than the force 
of gravity. Hence for this particular numerical example, they would go unnoticed by 
local gravity experiments [SI] , although, in general, couplings between the dark and the 
visible sectors can have observable consequences on structure formation [52"] . 

4. Conclusion 

We have entertained the possibility that the dark matter is composed by a Bose-Einstein 
condensate of light scalar particles coupled to the baryonic component of the Universe. 
The long-range forces among baryons caused by the pair exchange of dark matter 
particles in the presence of the condensate gives rise to a shift in the position of the 
acoustic peaks in CMB power spectrum above the first one. 

For a particular choice of the parameters in our model, a consistent cosmic evolution 
results that would show up as variations in the peak locations of the order of a few 
percent. Since our primary aim was to display the consequences of Bose-Einstein 
condensation of a relativistic scalar gas in a physical setting, we do not find it justified to 
exhaustively explore the full m<f>—g e ff parameter space. Other choices of the parameters 
could lead to a completely different, and even inconsistent, cosmic evolution. For 
example, the transition to a non-condensed phase could take place before the thermal 
decoupling, and we would be left with light scalar dark matter particles thinly populating 
the states above the ground state with no trace of the primordial condensate before 
structure formation occurs. Also, tests of the gravitational inverse-square law [SB] or 
Casimir force measurements j^lj, could be sensitive to parts of the — g e ff parameter 
space different from the ones studied. 

In the particular case considered, we find that the particles in the condensate will 
remain in the lowest energy state beyond the critical point. Once decoupled, the particles 
in the condensate are in practice isolated from the ones in the thermally excited states 
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and momentum redshifting can not cause any flow between this two populations. 
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